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Abstract 

We consider SUSY sine-Gordon theory in the framework of perturbed conformal 
field theory. Using an argument from Zamolodchikov, we obtain the vacuum structure 
and the kink adjacency diagram of the theory, which is cross-checked against the exact 
S matrix prediction, first-order perturbed conformal field theory (PCFT), the NLIE 
method and truncated conformal space approach. We provide evidence for consistency 
between the usual Lagrangian description and PCFT on the one hand, and between 
PCFT, NLIE and a massgap formula conjectured by Baseilhac and Fateev, on the 
other. In addition, we extend the NLIE description to all the vacua of the theory. 



1 Introduction 

Supersymmetric (SUSY) sine-Gordon theory (SSG) is of interest for several reasons. First, 
it is a two-dimensional integrable field theory and as such, there is a wealth of analytic and 
non-perturbative information available about its behaviour. Second, it describes supersym- 
metric solitons, which have been investigated recently in the context of supersymmetric 
gauge theories in four (and other) dimensions. It can be expected that SSG provides a 
useful laboratory for the analysis of solitons in supersymmetric theories, and their non- 
perturbative behaviour in general. Third, there has been a lot of activity recently especially 
concerning the situation when SSG has a nontrivial boundary condition imposed on a time- 
like boundary. 
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In spite of all the progress, much less is known about SSG than its non-supersymmetric 
counterpart, i.e. ordinary sine-Gordon theory. Although its S matrix has been conjectured 
several years ago [1], there is much less evidence for its correctness than in the non-SUSY 
case. In particular, the issue of quantum corrections to the soliton mass has been settled 
only fairly recently in [2], where it was proven that it satisfies a BPS-like property, even 
though the theory only has N = 1 supersymmetry. 

There is also a certain controversy in the literature concerning the vacuum and kink 
structure of the tricritical Ising model [3], the simplest nontrivial SUSY integrable field 
theory, which has a bearing on the S matrix of the SSG as well, since the S matrix of 
tricritical Ising model provides the description of the supersymmetric structure. 

In this paper we set out to clarify the vacuum and kink structure of the model (with- 
out boundaries), using information from several sources. First, in Section 3 we apply an 
argument going back to Zamolodchikov [4] which gives a kink structure compatible with 
the S matrix conjectured in [1]. In Section 4 SSG is considered in the perturbed conformal 
field theory framework, which is not entirely trivial, as a purely bosonic potential term has 
to be omitted from the Lagrangian (it is expected to be generated by radiative corrections 
in this description). To show that this description is correct we construct a spin-3 con- 
served charge and verify it against the classical limit. In addition, we analyze first-order 
perturbative corrections of the lowest-lying energy levels, and in Section 5 compare them 
to the results obtained from a conjectured NLIE, first derived in [5], which is extended 
to provide a finite volume description of all the vacuum levels. Finally, in Section 6 we 
use the truncated conformal space method to investigate the finite volume spectrum. Our 
conclusions are presented in Section 7. 

2 SSG theory 

2.1 Action and discrete symmetries 

The supersymmetric sine-Gordon (SSG) theory is defined by the action 

^4 SSG = J dtdx + flffd^ + m** cos ^$ + ^- cos p^j (2.1) 

where $ is a real scalar, ^ is a Majorana fermion field, m is a mass parameter and (5 is the 
coupling constant. The theory is invariant under an N — 1 supersymmetry algebra and 
has infinitely many commuting local conserved charges [6]. These charges survive at the 
quantum level and render the theory integrable, which makes it possible to describe the 
exact spectrum and the S matrix. We use the Weyl representation for the spinor field 

'■(-".;) ■ '-(:;) ■ '•-*'-(»':) 
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where ip± are real Weyl components with definite chirality. 

The supersymmetric theory has some discrete symmetries that play an important role 
in what follows. The field theory interaction is periodic 

An 

$^$ + n— , neZ (2.2) 

and even in the boson field 

$ -> . (2.3) 

These symmetries are also present in the non-supersymmetric version of the sine-Gordon 
model. It is also interesting to observe that the Lagrangian is invariant under a half-period 
shift 

2tt , x 

$ $ + -a ( 2 - 4 ) 

if at the same time one changes the relative sign of the fermion components, say 

^+^-^+ , i>-^i>- (2.5) 



or 



2.2 Spectrum and scattering amplitudes 

The spectrum consists of the soliton/antisoliton multiplet, realizing supersymmetry in a 
nonlocal way, and breathers that are bound states of a soliton with an antisoliton. 

The building blocks of supersymmetric factorized scattering theory were first described 
in [7], using an Ansatz in which the full scattering amplitude is a direct product of a part 
carrying the SUSY structures and a part describing all the rest of the dynamics. The full 
SSG S matrix was constructed in [1]. 

The supersymmetric solitons are described by RSOS kinks K € ah (0) of mass M and ra- 
pidity 0, where a, b take the values 0, \ and 1 with \a — b\ = 1/2, and describe the supersym- 
metric structure, while e = ± corresponds to topological charge ±1 (soliton/antisoliton). 
Multi-particle asymptotic states are built as follows 



Kl ai (*i) K\ a2 (02) ■ ■ ■ (^-0 K > 



(On) 



(2.6) 



where 9\ > 9 2 > . . . > 9 N for an in state and 9\ < 9 2 < . . . < 9 N for an out state. The 
two-particle scattering process 



Kl (*i) + Kc (02) - Kl &) + Kl (9,: 

has an amplitude of the form 



S. 



SUSY 



a d 
b c 



01 - 02 x S SG (0i - 02 , A)^ 



e l € 2 

ei«2 



(2.7) 
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i.e. the tensor structure of the scattering amplitude factorizes into a part describing the 
SUSY structure (which we call the SUSY factor) and another part corresponding to the 
topological charge (the bosonic factor). 

The bosonic factor coincides with the usual sine-Gordon S matrix, but the relation 
between the parameter A and (3 is different from the sine-Gordon case 

8tt 1 

The SUSY factor is identical to the S matrix of the tricritical Ising model perturbed by 
the primary field of dimension | [4]: 
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K{9) 



i fr £(fczi±A)£ftzA) 

v^r(Jb + i-£)r(Jb + £) 



= 2 e / 2 - C os(^ K{9) 
= 2 (--«)/^cos^ + ^ 
= 2^cos (1-|)^) 



(2.8) 



The supersymmetry factor has no poles in the physical strip and so the full supersymmetric 
amplitudes of the kinks have the same poles as in the ordinary sine-Gordon model. The 
corresponding bound states are supersymmetric breathers B n of mass 



Till 

m n = 2M sin — , n = 1, . . . , [A] 



(2.9) 



which transform in the ordinary boson-fermion doublet representation of the SUSY algebra. 
The bosonic component of B n has a nontrivial (— l) n parity under bosonic field reflection 



3 Kink and vacuum structure of SUSY models 

Here we use an argument by Zamolodchikov [4] to clarify the vacuum and kink structure 
of SUSY sine-Gordon theory. In its original form the argument is for the tricritical Ising 
model; we briefly recall it here. 
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Figure 3.1: The scalar potential in the Landau-Ginzburg description of the tricritical Isin£ 
model 
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Figure 3.2: The adjacency diagram of kinks in the tricritical Ising model 

3.1 Tricritical Ising model 



In the bosonic Landau-Ginzburg formulation, the basic field is the (bosonic) spin field a, 
and the action is of the form 

X - (d t a) 2 - l - (d x a) 2 - V») dxdt , 

where the potential has the form 

V(a) = X(a 2 -a 2 ) 2 a 2 . 

As can be seen from the form of the potential (Fig. 3.1), there are three vacua in the 
model, and the kinks connect them according to the adjacency diagram in Fig. 3.2, with 
a scattering amplitude given in (2.8). 

However, in the supersymmetric Landau-Ginzburg formalism the action (in component 
form) looks like 

\ {d&f ~ \ (<9*$) 2 + 4>+d^+ + i).d + i)_ + W'{§)i) + i)- + W{§) 2 j dxdt , 
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Figure 3.3: The scalar potential in the SUSY LG description of the tricritical Ising model 
where 

= n ($ 2 - A 2 ) . 
The scalar potential is of the form (Fig' 3.3) 

U($) =/ x 2 ($ 2 -A 2 ) 2 

and has two minima 

$ = ±A . 

How do we reconcile this picture with the bosonic one, in which there are three vacua? 
Zamolodchikov's argument goes as follows: at the minimum $ = A, the fermion mass is 
positive, there is only one vacuum, the Majorana fermion describes the high temperature 
phase of the Ising model and the vacuum expectation value of the spin field vanishes: 

(a) = . 

However, when $ = —A, we are in the low temperature phase of the Ising model. Therefore 
this vacuum is twofold degenerate, with 

(a) = ±a . 

The kinks go from the (a) = vacuum to the (a) = ±a ones. Therefore there are two 
pairs of kinks, and we get the same adjacency diagram (Fig. 3.2), with the (a) = (RSOS 
label |) vacuum in the middle, while the other two (ex) = ±<t (RSOS labels 0, 1) sit at 
the edges. 

It must be mentioned that Fendley has conjectured another bootstrap S matrix for the 
tricritical Ising model in [3], based on an assumption that the kink structure in the SUSY 
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Landau-Ginzburg formulation was different from that of the bosonic formulation. The 
bosonic and the super Landau-Ginzburg formulations are two different local descriptions of 
the model; their relation is very similar to the sine-Gordon/massive Thirring connection [8]. 
However, in the case of the sine-Gordon/massive Thirring relation the kink S matrix is 
independent of the local description taken, and Zamolodchikov's argument indicates that 
this is the case for the tricritical Ising model as well. As a result, in both the bosonic and the 
fermionic description the kink scattering is described by the same RSOS S matrix (2.8), 
contrary to Fendley's assumption. 1 The difference between the local algebras manifests 
itself in the different boundary conditions prescribed for the local fields, just as in the 
sine-Gordon/massive Thirring case (cf. also [9] for the situation in finite volume). 

3.2 N=l SUSY sine-Gordon theory 

Now let us apply these considerations to N — 1 SUSY sine-Gordon theory. Here the 
bosonic potential is 

- cos /3$ 

while the mass term of the fermion is of the form 

cos — . 

As a result, every even vacuum $ = 2n^ , n G Z is nondegenerate, while every odd 
vacuum $ = (2n + 1)^ , n e Z is doubly degenerate, with nonzero expectation value of 
the spin field associated to the Majorana fermion. We have a vacuum structure analogous 
to that of tricritical Ising model, which periodically repeats itself (similarly to the periodic 
vacuum structure of the sine-Gordon model). This results in the adjacency diagram shown 
in Fig. 3.4, which is exactly the one that describes the tensor product S matrix of kinks 
K±, a = 0,§,l. 

The discrete symmetries of the Lagrangian are manifest on the adjacency diagram, 
which is periodic in $ with period Att//3. The transformation (2.4, 2.5) is a translation 
with 2n/ (3, together with the analogue of Kramers- Wannier duality acting on the Majorana 
fermion, and exchanges the broken phase of the fermion (two vacua) and the symmetric 
phase (a single vacuum). 

Due to the periodicity of the vacuum structure and the field theoretic potential, in 
finite volume one can create twisted sectors depending on a i? angle according to [10], or 
one can define folded versions of the model by identifying vacua after a certain number of 
periods [11]. 

J TBA results also show that Fendley's S matrix yields an irrational value for the ultraviolet central 
charge, thus it cannot correspond to a perturbation of tricritical Ising model (M. Moriconi, private com- 
munication) . 
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Figure 3.4: A period of the adjacency diagram of kinks in supersymmetric sine-Gordon 
theory 

3.3 The IR description 

Now we give a qualitative description of the consequences of imposing periodic boundary 
conditions in the SSG theory in the IR domain, i.e in a large but finite volume L. It is 
based on the known features of the spectrum and the exact S matrix. 

The first consequence of imposing periodic boundary conditions is that the single kink 
states are excluded. Some of the two kink states are allowed; namely those K a bKb c where 
the a and c vacua are identical. As two kinks can fuse into breathers one may have single 
breather states in the sectors where the topological charge vanishes. It is easy to see, using 
the a — > 1 — a symmetry of the SUSY factor (2.8), that the allowed two kink states are 
divided into two sets, 

{K iKi , K lh K hl ) (Ki K oh , K hl K lk ) 

closed under scattering. The first set contains two vacua (the ones indexed by and 1), 
while the second one has a single vacuum (the one indexed by |). It is a natural expectation 
that in a finite volume these three vacua are no longer degenerate but their energies are 
exponentially close to each other. Above all three vacua we also expect the repetition of 
the spectrum of (single) breather states albeit with their energies split. One of the aims of 
this paper is to give a qualitative verification of this picture. 

The two sets of periodic two kink states have rather different supersymmetry properties. 
To establish this we use the kink-kink-breather fusion coefficients 

\K ab {6 + ic/2)K bc (o - <e/2)> = fL\m) + fL\m), 

(where \<p) and stand for the bosonic and fermionic breathers respectively) explicitly 
given in [12]. These explicit expressions reveal that the two kink states in the first set 
can fuse into bosonic breathers only, while the ones in the second set can fuse into both 
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bosonic and fermionic breathers. Furthermore, using the action of supersymmetry charges 
(also given in [12]), one can show that there is an unbroken N — 1 supersymmetry in the 
space of states of the second set. 



4 SSG as a perturbed conformal field theory 
4.1 The perturbed conformal field theory picture 

We consider SUSY sine-Gordon theory as a perturbation of the c = 3/2 conformal field 
theory of a free boson and a free fermion. The Hilbert space of the theory is taken to be 
a tensor product of the conformal free boson and free fermion, and we consider only local 
sectors, i.e. retain only states which have integer conformal spin. At generic values of the 
compactification radius this permits only two boundary conditions on the fermion field 

ty(x + L) = ±q?(x) 

which means either Neveu-Schwarz or Ramond boundary conditions on both components 
of the fermion fields (different boundary conditions on the components would result in 
sectors of noninteger spin) . 

We take a conformal boson field compactified on a circle of radius r, normalized with 
the action 

/ d 2 zd<\>d4> 



8vr _ 

which means that the canonical boson field of (2.1) can be written as $ = -j^<t>- The 
conformal boson has a (7(1) x U(l) symmetry generated by the modes of the currents 



J(z) = id(f) = a n z 11 1 and J(z) = id(p = a,, 



~—n—l 



The highest weight vectors \n, m) of the current algebra are created by the vertex operators 

In order to have a maximal local operator algebra we take both n and m to be integers. 
The Hamiltonian, on the cylinder of circumference L, in the sector \n, m) has the following 
form 

27r / /nv 



H boson L \{ r ) 



fc>0 / 



where the bosonic modes satisfy 

[a k ,ai\ = k5 k +i ; [a k ,a{\ = ; [a k ,ai] = kS k+ i 
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In the Neveu-Schwarz sector, the fermionic part of the Hilbert space can be generated 
by the negative modes of the Euclidean fermionic fields (corresponding to the Minkowski 
fields tp±) 

S+1/2&Z S+1/2SZ 

acting on the Neveu-Schwarz vacuum \NS) and locality requires the total fermion number 
to be even. The modes satisfy 

{b s ,b t } = 5 s+t ; {b s ,b t } = ; {b s ,b t } = 5 s+t . 
The fermionic Hamiltonian takes the form (using conformal normalization conventions) 

1 f 2,71 I °° i 

Hns = ^ / d 2 z tpdip = — b- s b s - — 

In the Ramond sector, the fermionic field has integer mode expansion and both the left 
and the right fermion number must be even. We denote the highest weight state by \R). 
The Hamiltonian in the Ramond sector is 

Writing cos f $ = |(V(i !0 ) +V(_i ; o)) gives the relation £ = between the compactification 
radius r and the coupling /3, and this leads to 

1 2 

P=T = — 7 • 

A r 2 — 1 




4.2 Conserved charges in the PCFT framework and the perturb- 
ing potential 

A convenient way to describe the SUSY sine-Gordon model is to consider it as an ap- 
propriate perturbation of the conformal field theory in its UV limit [13] by the operator 
\] 1 = V^cos ^. In [14] the SSG is described perturbation of the super Liouville theory 

2 

and it is shown that the purely bosonic piece of the potential U2 — cos (3$ vanishes in 
the renormalized theory. 

The problem is that assuming the kinetic terms of the boson and fermion fields in the 
Lagrangian provide the description of the c = | theory, the two parts of the interaction 

Ui = tpip cos @y and U 2 have different conformal dimensions: A l = A x = | + ^ and 
A 2 = A 2 = This, and the fact that they have different supersymmetry properties 
would complicate the PCFT description. Nevertheless - using the conserved quantities - 
we argue below that for the SSG it is enough to consider only U\ as a perturbation. 
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This perturbation obviously preserves a (1,1) supersymmetry generated by the super- 
currents G and G, since in the superfield formalism U\ can be written as G_i/ 2 G'_i/2 cos(|$) 
with cos(f $) being a Neveu-Schwarz superconformal primary field [15]. Thus the real ques- 
tion is whether U\ provides an integrable perturbation of the c = | theory. For this it is 
enough if U\ preserves a single higher spin conserved quantity. 

One can look for a conserved spin 3 quantity generated by the density T 4 , where 

t 4 = (d» 2 + A{d z <$>f + B(d z $) 2 d z ^ + Cd 2 z i,d z i, , 

with constants A, B and C. In the first order of the perturbation T 4 is conserved if in the 
operator product T±{z)Ui{w , w) the residue of the first order pole is a total derivative [13]. 
Having computed the residue one can eliminate the terms containing the derivatives of the 
ip field by using total derivatives; the remaining (non total derivative) terms then have the 
form 

(i^f$ + Bijjd 2 z <Z>d z <Z> + CV>(<9,$) 3 ) cos(^$) . 

Demanding the vanishing of A, B and C yields a linear inhomogeneous system of equations 
for A, B, C, that has a solution, which shows that Ui generates an integrable perturbation. 
Furthermore the (3 — > limits of A, B, C reproduce the classical expressions obtained in [16] 
using both Ui and Ui as perturbations together with the classical equations of motion. This 
provides the justification for only using U\ as the perturbation in the PCFT framework. 

Therefore in terms of the canonically normalized fields the Lagrangian of the SSG is 
written in the perturbed CFT framework as 

C = -d^d 11 ® + i^Yd^ + cos -$ . 

The coefficient of the perturbing potential is denoted here by \i to emphasize that its (f3 
dependent) dimension is different from that of the classical mass m. Identifying the bulk 
SUSY sine-Gordon model with the n = 2 case of [17] gives the following relation between 
the kink mass M and the \i parameter 

A* f l P 2 \ *,i (* P 2 , \ r(x) 

This relation - which connects an IR and a UV parameters and is called the massgap 
relation - plays an important role in writing the TCSA program devised to analyze the 
finite volume spectrum of SSG. Note that for (5 — > 0, the massgap relation correctly 
reproduces the classical kink (soliton) mass M — p. 

4.3 The PCFT analysis of energy levels in the UV 

The first PCFT correction to the UV energy levels can be expressed in terms of certain 
integrals over the complex plane of various 2p and 4p functions in the underlying conformal 
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field theory [18]. Indeed writing the (bare) Euclidean action of the perturbed CFT on the 
cylinder as 

A 9 = A CFT + 9 J d 2 Cx(0, 

where g has mass dimension y = 2 — 2A X and x( z ) is normalized by (x(z,z)x(0,0)) = 
|z|~ 4A x ; for small dimensionless volumes / = LM, the first correction to the energy of a 
state | a), with L \a) = A a \a) = L \a), can be written as 

-E a {l) = -(c - 24A a ) - c a 2 K 2 l 2y + o(l 4y ). 

TT 

Here c is the central charge of the UV CFT, it is assumed that the symmetry properties of 
the fields forbid any first order (o(/ y )) corrections and to derive this formula the massgap 
relation is written in the form g = kM v . We remark that the bulk energy constant in 
SSG vanishes on account of the supersymmetry. The coefficient of the leading correction 
at / — > can be written explicitly as 

eg = e^Tr) 1 " 2 " J 0(a| X (l,l)xM)|a)|conn. (4.2) 
plane 

Using the conformal normalized free fermions tjj, ip and the properly normalized V(±i,o) 
vertex operators to describe the perturbation in the SUSY sine-Gordon model one obtains 

X(z, z) = $(z)i/>(z)-j=(V( 1>0) (z, z) + V { _ lfi) (z, z)) 

and 

This perturbation is relevant if the compactification radius r > 1, and then |0), |1±), |2±) 
and |3±) constitute most of the lowest lying states in the various NS sectors. (Here \N±), 
N = 1,2,3 denote the states |JV±) = V (JV , )(0, 0)|0) ± V ( _ m (0,0)\0)) with |0) being 

the tensor product of the bosonic and NS vacua |0) = |0) <S> \NS)). For these states it is 
straightforward to determine the conformal correlation functions needed in eq.(4.2), and 
the integrals can be evaluated using the general expression 

/ d2z JT^ = *7 (" " f ) 7 (0 7(1 " ") ■ (4-4) 

(The integral is convergent at z — > if p > 0, at z ^ 1 if v < 1 and at z — > oo if 
2v> — p > 0; we keep the conditions 2v> > p > but continue in v above v — 1 using the 
r.h.s. if necessary) . The coefficients of the leading PCFT corrections for these seven states 
are summarized in Table 4.1, where the coefficient a is 

3, . 2 1 / 1 \ 
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State 


A a 




|0> 





«7 2 (| + 7 (-£) 




1±) 


1 

2r 2 


«7(i - W-^)7(i + ± HI - ^)7(^)] 




2±> 


2 

r 2 


«7(|-^)7(-^)7(| + ^) 




3±) 


9 
Ir 2 " 


«7(i-^)7(-^)7(i + ^) 



Table 4.1: The first PCFT corrections for some states 

Note that in all these cases it is necessary to make the analytical continuation in v 
to give meaning to the otherwise divergent integrals. This divergence is expected as the 
conformal dimension of the perturbing operator A = | + ^2 is greater than 1/2; thus the 
singularity coming at z — > 1 from the x(^)x( z ) OPE is not integrable no matter what the 
states I a) are. 

In the Ramond sector we expect that the ground state scaling function vanishes identi- 
cally as a result of unbroken supersymmetry. The question is whether the (leading) PCFT 
corrections are consistent with this expectation. The conformal contribution to the scaling 
function vanishes since 

c - 24A R =-- 24— = . 
H 2 16 

The next contribution is proportional to the second order coefficient 

=6(27r) 1 - a » J ^(R\ X (l,l)x(z,z)\R)\ 



where 



conn 

\z\» 
plane 

\R) = lim a(z,z)\NS) 



denotes the Ramond ground state. The integrand defining c 2 is formally positive, but 
the integral is divergent and needs regularization. Separating the integrand into a sum of 
terms, after an appropriate partial integration it can be converted into a form in which each 
term can be regularized using (4.4). The sum of regularized terms turns out to vanish. 
In this sense PCFT is consistent with a vanishing ground state scaling function in the 
Ramond sector. 



5 NLIE 

To probe the behaviour of the scaling functions between the ultraviolet and infrared limits 
one can use the thermodynamic Bethe Ansatz [19,20]. Fendley and Intriligator's TBA 
gluing idea [21] leads to a set of TBA equations for all choices of the SSG coupling [22-24], 
but in practice the equations take a simple form only at certain (3. Instead one can treat all 



13 



values of the coupling on an equal footing by using an alternative type of nonlinear integral 
equation, originally developed for the ground state of the sine-Gordon model in [25] (and 
independently in a related context in [26,27]), and which usually goes by the acronym 
NLIE. 

The equations proposed in [5] for the SSG model consist of a TBA-like function coupled 
to a nonlinear integral equation. The two pieces reflect the factorized nature of the S matrix, 
with the TBA part associated to the SUSY factor and the NLIE part to the bosonic factor. 
The original equations of [5] generate the lowest ground state in the Neveu-Schwarz sector, 
but a simple sign change and appropriate choice of parameters provides access to the other 
two ground states. The equations are 

/oo 
d9' X (0-9'+f )ln(l + (-l) 5 y 2 (#')) 
■oo 

+ / d9'<p(9-9')\n(l + yi (9'))- [ d9'<p(9-9')\n(l + y^(9')) ■ 
JCi Jc 2 

\ny 2 (9) = [ d9'x(9-9'-f)Hl + yi (9'))- [ d& *(0-0'-f ) ln(l + y?{9')) . 

The integrations contours C\ and C 2 run from —00 to 00 just below and above the real 
axis respectively. The kernel ip{6) is proportional to the logarithmic derivative of the 
soliton-soliton scattering amplitude of the sine-Gordon model 



dk ike sinh(p-l) 



2 



2sinh^coshf ' 
and x(Q) i s related to the SUSY factor of the S matrix 

X(9) = l/(27rcoshfl). 
The exact ground state energy depends only on the NLIE-like function via y\{9) 

E(l) = -7^ ^ d9 sinh^ \n(l+ yi (9)) -J d9 sinh 9 ln(l + yf 1 ^))^ . 

By setting the parameters uj and 5 appropriately these equations provide access to all 
three vacua, though not for all possible choices of (f3, 1) for each state, as we shall comment 
on further below. The appropriate settings are shown in table 5.1. 

Typically we can extract a number of results analytically from nonlinear integral equa- 
tions, the simplest of which is the ultraviolet value of the ground state energy. In the 
Neveu-Schwarz sector (5 = 0) the scaling function behaves as : 

61 



61 , . /3 12u; 2 \ 

= -(2-—) 
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State 


a; 


5 


|0) <g> \NS) 








|l+> 


1 





|0) <g> 


1 

2 


1 



Table 5.1: Choice of NLIE parameters to obtain each of the three ground states. 



while in the Ramond sector (5 = 1) we have 



61 „, . / 12(w- 1/2)' 
—E(0) = 1 



7r \ r 2 

Note that tuning cu according to Table 5.1 yields the correct UV behaviour for each ground 
state: 

-E(0) = - (c - 24A a ) . 

71 

Furthermore the conformal dimension of the perturbing operator exactly matches that of 
U\. In the infrared limit (/ — >• oo) we deduce the scaling functions behaviour to be [5] 



/d9 
— cosh 6 e 
2tt 



Ml cosh 8 



(5.1) 



This result can be given a simple intuitive interpretation. First, let us note that the 
RSOS structure describing the SUSY degrees of freedom of solitons has three nodes 0, 1/2, 
1, which correspond to truncating the quantum group representation theory U q (sl(2)) at 
g 4 = 1 to the tensor product rule 

0®a = a , a = 0,-,l 

\®\ = §©\ , 101 = . 

There is a notion of statistical dimension for these representations d a (a = 0, 1/2, 1), which 
satisfies the ordinary rules under the tensor product 

1 

«o«a — d a , a = 0, -, 1 
d\ — d + di , d\ — d , 

2 

and the solution to these conditions is 

do — d\ — 1 , di = V2 . 

The kinks correspond to the 1/2 representation, and the above results mean that they have 
a statistical dimension \/2. It might seem strange since it is hard to interpret this as the 
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r 


b 2 


c%k 2 


>/l5/3 


-0.4305230886667 


-0.4305230886637 


v/21/3 


-0.7099942272385 


-0.7099942272382 


2 


-0.5922281367467 


-0.5922281367458 


3 


-0.2517882823698 


-0.2517882823567 


3.05 


-0.2422117898616 


-0.2422117898994 



Table 5.2: Comparison of PCFT correction with NLIE for the first Neveu-Schwarz ground 
state 

multiplet length. However, a simple calculation shows that the number of 2n-kink states 
(neglecting for the moment the topological charge) with periodic boundary conditions grows 
as 2 n (there are no odd-kink states on a circle) so the kinks indeed behave as a multiplet of 
length a/2- The reason for this can be found in the nontrivial vacuum adjacency conditions 
imposed on multi-kink states, which do not allow for all possible sequences to be realized. 

Armed with this, it is easy to interpret the result (5.1). For u = it is exactly the 
leading correction from 2\/2 particles in finite volume (or equivalently finite temperature) 
to the free energy. The additional factor 2 comes from the fact that the topological charge 
introduces a further doublet structure of the kinks. The uj dependence is easy to understand 
using the same sort of instanton argument as in the case of the ordinary sine-Gordon 
model [11]. 2 

A final check can be made by comparing NLIE data with the above perturbed conformal 
theory predictions. Turning to the first Neveu-Schwarz ground state, we extracted the 
coefficient b 2 from NLIE data fitted to a suitably truncated series of the form 

«/ 

-E a (l) = -(c - 24A a ) -VW 2 " 9 , 

ti 

then made a comparison with the coefficient of the leading PCFT correction c^k 2 . The 
results displayed in Table 5.2 show excellent agreement. Numerical iteration of the NLIEs 
in the Ramond sector confirms the Ramond ground state energy is identically zero for all 
choices of This provides a nontrivial check of the equations since the NLIE functions 
yi{0) and y 2 (9) are not identically zero. 

It turns out to be harder to check the NLIEs for the second Neveu-Schwarz ground 
state. For small values of r and large values of the cylinder size / the equations work well. 
This allows a comparison with PCFT, but does not allow one to obtain the value of c%k 2 to 

2 We remark that a similar explanation can be made for the IR behaviour of the vacuum scaling function 
predicted by the NLIE, for all the fractional supersymmetric sine-Gordon models, c.f. formula (29) in 
[5]. For the case L of that paper, the statistical dimension of the kink representation 1/2 results in 
2cos(tt/(L + 2)). SSG corresponds to L = 2, and ordinary sine-Gordon to L = 1 (i.e. the kinks have 
dimension 1 apart from the double degeneracy corresponding to their topological charge) . 
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Figure 5.1: Comparing PCFT with NLIE for the second Neveu-Schwarz ground state 

as many decimal places as for the lower NS ground state. Instead, in Fig. 5.1 we compare 
the PCFT result at compactification radius r 2 = 11 (equivalently p = 1/5) 

c-24A| 1+> + c[ 1+) k 2 l 2y , 

using a solid line to the NLIE data —6lE a (l)/ir (indicated by the symbol o). The agreement 
is very good. 

So what lies behind the difficulties at some couplings and/or cylinder size? For larger 
values of r or smaller / the equations fail because one or more zeros of 1 + y\{9) attempt to 
cross one of the integration contours, resulting in a singularity in ln(l + y 1 (9)). (Once this 
has happened a zero of 1 + 1/2(0) may also attempt to cross a contour.) This behaviour is 
not unexpected: by sending the cylinder size / into the complex plane the authors of [28,29] 
were able to analytically continue a set of TBA equations describing the ground state energy 
of a particular model to a new set, which gave access to the first excited state. Here, we 
have an additional parameter 00 which we analytically continue from to 1 along the real 
axis to obtain the second NS ground state from the first. Singularities in the logarithmic 
term of NLIEs of the type discussed here also arise via a second mechanism, which usually 
occurs as the (real) cylinder size / is decreased and is not explicitly connected to an excited 
state. This so-called 'specials' problem was first described in the context of the ordinary 
sine-Gordon NLIE when the equation was tuned to study the various states of the minimal 
models perturbed by 0i 3 [30]. Since we are able to make a favourable comparison of NLIE 
data and PCFT results for some couplings r we leave the resolution of these problems for 
future investigation. This will also provide clues as how to modify the NLIEs to address 
the excited states, as has been done for the sine-Gordon NLIE in [9,31-34]. 
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6 TCSA 



6.1 TCSA for the supersymmetric sine-Gordon model 

The supersymmetric sine-Gordon model can be viewed as a relevant perturbation of the 
c = | free conformal field theory, consisting of a free fermion and a free boson. Following the 
idea of the truncated conformal space approach (TCSA) [35], the Hilbert space is truncated 
at a given conformal energy E cu ^ , and so the Hamiltonian can be diagonalized numerically 
to obtain an approximate finite volume spectrum. Since the perturbing operator has scaling 
dimension larger than 1/2 we are faced with divergences in the TCSA data and hope to 
obtain relative energy levels only. Furthermore, the Hilbert space contains many more 
vectors up to a given energy cut than in a theory with a single scalar field only, so only 
qualitative results are expected. 

The matrix elements of the perturbing operator between any two conformal states \a) 
and \b) with conformal weights (A a , A a ) and (A&, A&) can be brought to the following form 

(a\gJ L dx X (tA\b) = ^ (jfj <a|iV(l)^(l) (V (1>0 )(1, 1) + V ( -i >0 )(l, 1)) \b) S Aa _ AaAb _ Ab 

using a conformal mapping from the cylinder to the plane and translational invariance (the 
S expresses the conservation of momentum/conformal spin). Recall that A = | + y?. 

We use relations (4.1,4.3) to write the PCFT coupling g in terms of the kink mass 
M and to convert the Hamiltonian into dimensionless form (energy measured in units of 
M), parametrized by the dimensionless volume parameter / = ML. The Hilbert space 
can be decomposed into sectors of given conformal spin and given bosonic winding number 
(topological charge). 

Using the mode expansions and U(l) charge conservation, the matrix elements of the 
perturbing operator can be calculated as the product of separate bosonic and fermionic 
contributions. In the Ramond sector, matrix elements including the fermion zero mode 
can be computed using the relation 

(R\ib b \R) = X - . 

Since the interaction is even in the bosonic field, the perturbing operator has the following 
discrete bosonic Z 2 symmetry 

| n, to) <-> | — n, —to) ; a n <-> — a n ; a n <-> — a n . 

The discrete symmetry (2.4,2.5) means that the NS sector has an additional Z 2 symmetry 
acting as 

| n, to) <-> (-l) n |n,m) ; b s <-> -b s ,b s ^b s . 

There is no such symmetry in the Ramond sector as a result of the presence of the fermionic 
zero mode. 
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|i-)r 


|2+) 


13+) 






|i+)r 





Table 6.1: UV classification of the lowest lying states in the various NS and R sectors for 
r > 3. For r < 3 the order of the states is somewhat different (|3+) is the third state in 
its sector), but they can still be identified unambiguously. 

The Hilbert space (see Section 4.1) can be further decomposed with respect to these 
symmetries. In table 6.1 we list the lowest energy states of each sector for zero spin and zero 
topological charge. The first sign refers to the NS parity, while the second to the bosonic 
one. In parenthesis we indicate the symbol we used in the TCSA data for the various 
sectors and |7V±) R = ^ (V^O, 0)|0) R ± V^o^O, 0)|0) R ), where |0) R = |6> ® \R). 

6.2 Checking TCSA against PCFT 

The aim of this investigation is to confirm the TCSA program by comparing the small vol- 
ume (UV) data for the low lying energy levels to the predictions of conformal perturbation 
theory (PCFT). Since the TCSA data are given in terms of the (IR) kink mass while the 
PCFT predictions in terms of the (UV) parameter /z, to make the formal comparison one 
has to use the massgap relation. Nevertheless as the TCSA program is written using the 
same massgap relation the success of the comparison says nothing about the correctness 
of this relation but may confirm the TCSA program. 

Since the scaling dimension of the perturbing operator is larger than 1/2 both the 
conformal perturbation theory and TCSA are plagued by divergences. As discussed in 
section 3.3, in PCFT these divergences are regularized by analytical continuation. The 
introduction of a finite E cu ^ in TCSA also regularizes these divergences albeit in a different 
way. Therefore the best one can do to compare TCSA and the PCFT predictions is to 
consider the differences between the various energy levels rather than the levels themselves: 

6/ 

-(E b (l) - E a (l)) - 24(A 6 - A a ) = Ac ba K 2 l 2y + o(Z 4 ^), Ac ba = c\-c\. 

7T 

This has two advantages: on the one hand the TCSA energy differences depend less sensi- 
tively on £" cu t than the individual levels, and on the other the integrals defining Ac ba are 
convergent thus they may be computed as the appropriate differences of the data in the 
table 4.1. 

Using the identifications between the low lying TCSA lines and the \N±) states in table 
6.1 we compare the TCSA energy differences and the PCFT predictions in Fig. 6.1. In 
diagram (a), the three continuous lines depict the PCFT predictions with Ac' 1 ^' ' 1 Ac' 1+ ^ 
and Ac' 2+ ^ respectively, while - in accordance with table 6.1 - the symbols *, o, □ and 
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x denote the TCSA data for 
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respectively. Note that the data is consistent with the |2±) states being degenerate in the 
leading order of PCFT. 

In diagram (b), the first three lines and the first four sets of TCSA data are just like 
on (a), while the fourth continuous line depicts the leading PCFT prediction for Ad 3+ ^ 1+ ^ 
and the second line of points marked by o are the TCSA data 

(£-+(,)-£-(,)) 

The agreement between the data and the leading predictions is excellent, though the de- 
viations between them indicate that in some cases the higher order PCFT corrections are 
not negligible even in this / region. 



6.3 IR spectrum from TCSA 

Since the scaling dimension of the perturbing operator is larger than 1/2 we have diver- 
gences in the TCSA data. This can be seen by increasing the cutoff: the energy levels 
become more and more negative. The energy differences, however, converge but unfortu- 
nately very slowly. The first two lines in the (++), ( — h) NS and in the (+) R sectors, 

and the first (lowest) lines in the (H — ), ( ) NS and (— ) R sectors are shown in Fig. 6.2, 

where we plot these lines relative to the ground state of the (++) NS sector. 

The various points are marked according to table 6.1. The first row in the table corre- 
sponds to the UV limit of the three vacua as shown on the data. The next row, (the next 
three curves) corresponds to the first breathers excited above each vacua, while the last 
row, (the last three curves) are related to the second breathers. The bosonic breathers are 
identified on the basis of their behaviour under the discrete symmetries. The compactifica- 
tion radius is r = 3.05 and the energy cut is 14, which corresponds to approximately 13000 
states in the Ramond sector and 9000 states in the NS sectors, respectively. Clearly the 
curves are in a qualitative agreement with the spectrum, but the particular mass values 
are wrong. Indeed it can be seen that the TCSA converges very slowly to the right values: 
by changing the truncation level from 8 to 14 in steps of 2, the energy differences change 
almost the same amount in each step. However, the qualitative features of the spectrum 
(i.e. the degeneracy and gap patterns) already appear for very small number of states and 
are stable against increasing the truncation level. 
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Figure 6.1: Comparison of TCSA and PCFT energies relative to the NS (++) vacuum 
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Figure 6.2: Lowest lying states in the TCSA spectrum 



Even though the mass values converge very slowly, the case for the mass ratios is much 
better. Calculating particular mass ratios from the data, the predictions of S matrix theory 
(2.9) can be reproduced with a few percent accuracy, and are much more stable against 
changing the cut. 

Furthermore, the TCSA data are also qualitatively consistent with the dependence of 
the spectrum on the coupling, i.e. by changing the compactification radius r the number 
and relative position of breather line triplets follows the predicted spectrum. 



7 Conclusions 

In this work we investigated supersymmetric sine-Gordon theory (SSG), especially the 
structure of the vacua and the kinks connecting them. Using an argument due to Zamolod- 
chikov we clarified the vacuum and kink structure and showed that it corresponds exactly 
to the exact S matrix conjectured in [1]. Then we formulated the theory in the perturbed 
conformal field theory (PCFT) framework, which shows a striking difference to the usual 
Lagrangian description, namely the omission of a purely scalar potential term. The first 
evidence we gave for the correctness of this description came from the construction of a 
spin-3 conserved charge, which is necessary for integrability. Considering the classical limit, 
we showed that the expression for this charge obtained in the PCFT framework reduces to 
the result known from the classical Lagrangian approach, evidence that the two approaches 
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indeed describe the same model. 

Using PCFT, we obtained the leading behaviour of some energy levels in finite volume, 
for the limit of small volume ( the ultraviolet regime). In particular, we concluded that 
the energy of the Ramond ground state vanishes in this limit, which is consistent with the 
fact that it is expected to vanish exactly due to unbroken SUSY in finite volume in this 
sector. 

The PCFT results were then compared to results obtained from an NLIE equation for 
the lowest lying Neveu-Schwarz state, which is the true ground state in finite volume. This 
comparison involves the mass gap relation between the PCFT coupling (an UV parameter) 
and the soliton mass (an IR quantity) conjectured earlier in [17]. The comparison showed 
excellent agreement for several different values of the dimensionless coupling parameter (5 
(equivalently r) of the model, which is a strong evidence for the correctness of the mass 
gap formula on one hand, and for the consistency between the PCFT and NLIE approach 
on the other. 

We have also examined the large volume (infrared) limit of the NLIE and found a 
perfect agreement with considerations based on statistics of kinks and instanton calculus 
in finite volume. 

The NLIE proposed in [5] was formulated for the true ground state of SSG, which is the 
lowest lying state in the Neveu-Schwarz sector. Matching conformal dimensions computed 
from the NLIE with those predicted from conformal field theory, we proposed equations 
describing the other two vacua, one in the Ramond and the other in the Neveu-Schwarz 
sector. Numerical iteration of the NLIE for the Ramond ground state gave a vanishing 
energy for all volumes, as expected, providing a very strong check since the solution of 
the equation itself is still nontrivial in this case. For the second Neveu-Schwarz vacuum, 
convergence properties of the iteration are not as good as for the other two cases, but 
the UV behaviour matches the PCFT prediction nicely, giving another confirmation of the 
consistency between PCFT, NLIE and the massgap formula. 

Finally, we performed an analysis of the finite volume problem using truncated confor- 
mal space approach (TCSA). We checked the validity of the approach in the UV regime 
by a comparison with PCFT. While the numerical convergence of the TCSA is not very 
good (partially due to UV divergences, but also to a rapidly growing number of states as 
the energy cut is raised), it does confirm the qualitative picture of the spectrum, both the 
degeneracies and the behaviour of the breather levels are in accordance with the picture 
presented before. 

There are many open issues. It would be very interesting to give a complete charac- 
terization of the SSG finite volume spectrum from the NLIE, like it has been done in the 
non-SUSY case [9,25,26,30-34]. It is not at all obvious how to achieve that, however, 
mainly because the NLIE for the SSG theory is of a mixed type - part of the system resem- 
bles the Destri-de Vega equation, while the other part looks like a thermodynamic Bethe 
Ansatz equation. It is likely that a reformulation of the NLIE is necessary. 

Another interesting problem would be to extend the analysis of the spectrum to the 
boundary case, along the lines of non-SUSY sine-Gordon theory [36]. This was part of the 
motivation to introduce TCSA for SSG, however, the slow convergence of TCSA prevents 
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us from obtaining reliable quantitative results, which would be necessary to check the 
bootstrap spectrum proposed in [37]. It remains to be seen whether this can be overcome 
in some way. Extending the NLIE to the boundary situation could be another way to 
obtain interesting results and check the bootstrap predictions. 
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